*****************************
** Main manuscript results **
*****************************


** this script reproduces the main model **
** results given in Tables 2 and 3 **
** (be sure to set your directory as appropriate) **

	use "data.dta", clear
	
* Table 2
	reg ihsPink lnseats  if democ == 1, robust cluster(cp)
	estimates store m1
	reg ihsPink lnseats gdp_growth gdp_capita eci_value female_labor_ilo								if democ == 1, robust cluster(cp)
	estimates store m2
	reg ihsPink lnseats gdp_growth gdp_capita eci_value female_labor_ilo Presidential left_exec pr 		if democ == 1, robust cluster(cp)
	estimates store m3
	reg ihsPink lnseats gdp_growth gdp_capita eci_value female_labor_ilo sd_ahs  						if democ == 1, robust cluster(cp)
	estimates store m4
	reg ihsPink lnseats gdp_growth gdp_capita eci_value female_labor_ilo pink_elasticity i.pair 		if democ == 1, robust cluster(cp)
	estimates store m5
	reg ihsPink lnseats gdp_growth gdp_capita eci_value female_labor_ilo oil_rents						if democ == 1, robust cluster(cp)
	estimates store m6

	estout m1 m2 m3 m4 m5 m6 , order(lnseats ) cells(b(star fmt(a2)) ///
		se(par fmt(%4.3f))) starlevels( * .1 ** .05 *** .01) stats(N , fmt(0) labels("Number Obs." ))  ///
		nolz style(tex) drop(*.pair) ///
		varlabels(lnseats "Log seat share women" gdp_growth "GDP growth" gdp_capita "GDP per capita" ///
		sd_ahs "Tariff dispersion" _cons "Constant" ///
		oil_rents "Oil wealth" female_labor_ilo "Female labor participation" ///
		pink_elasticity "Difference elasticity" left_exec "Left Chief Excutive" ///
		eci_value "Economic complexity" pr "Proportional representation"  ///
		L.pinkTax_ahs "Lagged DV" Presidential "Presidential system"  ///
		lnimport_men "Log imports men's item" lnimport_women "Log imports women's item" )

* Table 3:
	reg ihsPink c.lnseats##i.nondemoc gdp_growth gdp_capita eci_value female_labor_ilo, 		robust cluster(cp)
	estimates store m1
	reg ihsPink lnseats gdp_growth gdp_capita eci_value female_labor_ilo i.year					if democ == 1, robust cluster(cp)
	estimates store m2
	reg ihsPink lnseats gdp_growth gdp_capita eci_value female_labor_ilo i.iso3n				if democ == 1, robust cluster(cp)
	estimates store m3
	reg ihsPink ihsPinkLag lnseats gdp_growth gdp_capita eci_value female_labor_ilo i.year 		if democ == 1, robust cluster(cp)
	estimates store m4
	reg ihsPink lnEnforceCount gdp_growth gdp_capita eci_value female_labor_ilo i.iso3n i.year	if demSample == 1, robust cluster(cp)
	estimates store m5
	reg ihsPink lnEnforceCount gdp_growth gdp_capita eci_value female_labor_ilo i.iso3n i.year	if demSample == 1 & ever_adopt == 1, robust cluster(cp)
	estimates store m6
  
	estout m1 m2 m3 m4 m5 m6 , order(lnseats c.lnseats lnEnforceCount) cells(b(star fmt(a2)) ///
		se(par fmt(%4.3f))) starlevels( * .1 ** .05 *** .01) stats(N , fmt(0) labels("Number Obs." ))  ///
		nolz style(tex) drop(*.iso3n *.year) ///
		varlabels(lnseats "Log seat share women" gdp_growth "GDP growth" gdp_capita "GDP per capita" ///
		lnEnforceCount "Years of quota Implementation" _cons "Constant" ///
		oil_rents "Oil wealth" female_labor_ilo "Female labor participation" ///
		eci_value "Economic complexity" pr "Proportional representation"  ///
		ihsPinkLag "Lagged DV" Presidential "Presidential system"  ///
		lnimport_men "Log imports men's item" lnimport_women "Log imports women's item")
